Measuring every particle's size from three-dimensional imaging experiments 
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Often experimentalists study particulate samples that are nominally monodisperse. In reality 
many samples have a polydispersity of 4-10%. At the level of an individual particle, the consequences 
of this polydispersity are unknown as it is difficult to measure an individual particle size from 
images of a dense sample. We propose a method to estimate individual particle radii from three- 
dimensional data of the particle positions. We validate our method by numerical simulations of four 
major systems: random close packing, colloidal gels, nominally monodisperse dense samples, and 
nominally binary dense samples. We then apply our method to experimental data from moderately 
concentrated colloidal suspensions observed with confocal microscopy. We demonstrate that we can 
recover the full particle size distribution in situ. Lastly, we use our method to study the relationship 
between homogeneous colloidal crystal nucleation and particle sizes. We show that nucleation occurs 
in regions that are more monodisperse than average. 
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A wide variety of techniques exist for three- 
dimensional imaging of collections of particles [l], 
These types of samples include granular materials, soil 
mechanics, and colloidal suspensions. Our particular in- 
terest is in colloidal suspensions; these have been suc- 
cessfully used as model systems for understanding phase 
transitions for several decades [H, H| , and moreover are 
interesting in their own right due to industrial relevance 
Q. Confocal microscopy can be used to take three- 
dimensional images of fluorescent colloidal particles deep 
within a sample [H, 01 ■ When coupled with particle 
tracking techniques, the motion of thousands of individ- 
ual colloidal particles can be followed over long periods of 
time |§4ll|. This technique has been used to investigate 
the colloidal glass transition [H[, crystallization 

EMI, colloidal gels EMI, capillary waves El [IT 
sedimentation [R 
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and a variety of other questions 
(see ref. El f° r a review). One advantage of confocal mi- 
croscopy of colloids is that the data obtained are similar 
to what is found using simulations, which also provide 
the data of particle positions over long periods of time. 

However, experimental samples are always polydis- 
perse: even for a nominally single-component sample, the 
particles have a variety of sizes [22| . This is quantified by 
the polydispersity p, defined as the standard deviation of 
particle radii divided by the mean radius. For many sam- 
ples, p <~ 0.04 — 0.10 f22|. From numerical simulations, we 
know that the effects of the particle size distribution are 
not negligible. For example, crystal nucleation is difficult 
or impossible for more polydisperse samples [23], [24| . The 
crystal-liquid phase boundary depends on the polydis- 
persity [251 ] . The sensitivity to volume fraction near the 
glass transition depends on the composition in nontrivial 



ways [26|, l27[ . Experimentally, the influence of polydis- 
persity on colloidal crystallization has been demonstrated 
[28l 129( 1 , and there is also some understanding of how the 
particle size distribution influences the rheological behav- 
ior of a colloidal sample [30]. However, these are limited 
to studies of the spatially averaged properties of the sam- 
ple. Microscopy is useful for local properties, but particle 
size fluctuations of 0.04 — 0.10 are not easily detectable. 
It would be desirable to know particle sizes for more di- 
rect comparison with simulations. Furthermore, in some 
cases, neglecting these sizes in an experiment can lead 
to wrong conclusions. One example is that the pair cor- 
relation function g(r) can show a qualitatively incorrect 
dependence on control parameters if the particle sizes are 
treated as all identical [3l|] . A second example is that the 
apparent compressibility of a random close packed sam- 
ple depends qualitatively on whether individual particle 
sizes are taken into account I32H34I . 
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In this work, we introduce a general method for using 
3D data to determine the size of individual particles in 
any moderately concentrated sample, in general with vol- 
ume fractions <f> > 0.4. Wc use simulation data to verify 
that our method works well in a variety of sample types. 
We then demonstrate the utility of our method using 
previously published experimental data from confocal mi- 
croscopy of colloids. In particular, we show that colloidal 
crystal nucleation is sensitive to the local polydispersity: 
nucleation happens in locally monodisperse regions. Our 
method is not limited to confocal microscopy and col- 
loidal samples, but rather works with any data of the 3D 
positions of a collection of particles. 

Due to diffraction limits, it is difficult to directly de- 
termine the radii of individual particles from microscopy 
images to better than ±0.1 /im [35). Defining the edge 
is somewhat arbitrary and varies depending on particle 
properties and the details of the microscope illumination. 
Other 3D imaging techniques have similar issues 0. 



2 



In contrast, it is much easier to calculate the mean ra- 
dius a of particles with a variety of techniques [22| . Like- 
wise, from the centers of particles, the separations be- 
tween neighboring particles can be easily calculated. 
Our estimation method for particle sizes uses only a and 
Tij. The key idea of our method is that a large particle 
will be slightly farther from its neighbors and thus have 
larger values for its r^ , and likewise a smaller particle 
will have smaller values of ry. 

To start, we relate the pairwise separations as 

r. tj (t) = di + a.j + Sij(t), (1) 

where particle j is a nearest neighbor particle of parti- 
cle i, Tij is the measured distance between i and j, a% 
and dj are their radii, and Sij(t) is a surface-to-surface 
distance between their particles. We typically consider 
5 — 7 nearest neighbor particles (the closest neighbors); 
this choice is justified below. Often these data come 
from particle tracking [ll], HH an d so r,j (t) and Sij (t) 
depend on time t. Next we take an average of with 
respect to the nearest neighbor particles j, and then 
(rij(t))j = di + (a,j)j + (Sij{t))j, where ()j means an 
average over particle j. Thus, we obtain 

a i = (nj{t))j - (aj}j - (5ij(t)}j. (2) 

This is exact, but the quantities 5ij(t) are unknown. 
We estimate this by replacing <%(t) with its time- and 
particle-averaged value, the mean gap distance 5 = 
( r ij(t))i,j.t ~ 2a, where the average is over all particle 
pairs and all times. Our algorithm is then: 

of (t) = a, (3) 
a { ? ] (t) = (rvWh-iaf-VWh-l (4) 

where the superscripts denote iteration. The more we 
iterate Eq. |4l the more information we obtain from par- 
ticles far away from a given particle. In fact, a\ (i) is 

unchanged for n > 10 since aj 10 ^ includes the informa- 
tion from several thousand particles, thus we fix n = 10 
for the number of iterations in this paper. Of course, 
the particle radius does not depend on time, so after the 
10th iteration, we time-average a\ (t) to obtain the es- 
timated particle radius a^ 10 \ Time-averaging after each 
iteration of Eq. @]does not change the results. 

There are several sources of uncertainty in this esti- 
mation. First, there is the uncertainty of each particle 
position. Typically this is about 5-8% of the mean radius, 
leading to a 8-10% uncertainty of [§, [ll], [3(|. How- 
ever, these errors are nearly time-independent, so those 
errors are greatly diminished by time averaging. Second, 
our approximation for 6 is weaker in the case that the 
distribution of Sij(t) is broad. This in part depends on 
how many nearest neighbor particles are chosen: more 
neighbors results in a broader distribution, whereas too 
few neighbors means that the average (rij(t)) in Eq. @]is 
poor. Below, we use simulation data to determine that 
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FIG. 1: (Color online) Dependence of estimation uncertainty 
Aa on parameters. (A) The uncertainty Aa as a function of 
the chosen number of neighbors Z used for the averaging. The 
circles, triangles and squares correspond to Aa for volume 
fractions 4> = 0.51, 0.56, and 0.64 (RCP), respectively. In 
each case, the sample polydispersity is p — 0.07. (B) Aa as 
a function of the bulk polydispersity p. The circles and the 
triangles correspond to Aa at RCP and cf> = 0.51, respectively. 
The solid lines are the fitting lines for Aa. 



Z = 5 — 7 nearest neighbors is an optimal choice. Third, 
independent of a given choice of Z, some particles will 
simply be farther from their neighbors, and some will be 
closer. In a dense suspension, for example, this relates 
to the size of the "cage" formed by the nearest neighbor 
particles [l(J. Again, time averaging helps. If particles 
can rearrange and find new neighbors, then S becomes a 
better approximation for {Sij(t))t- In dense colloidal sus- 
pensions with volume fractions <fi > 0.5, rearrangements 
become infrequent and so longer time averages are de- 
sired [8T4l0|. In summary, the greatest strength of our 
algorithm is time-averaging, and past that, a sensible 
choice for the number of nearest neighbors Z is useful. 
Our tests show that time averaging over ~ 20 different 
times is sufficient for reasonable results. 

To verify our radius estimation method, we simulate a 
variety of systems and compare the estimated radius of 
each particle with its true radius. The error is given by 
5a,i = a| 10 ' 1 — di, where a| 10 ' is the estimated value and 
di is the true value. Aa = {5di)/d is the mean fractional 
error in the estimated particle radius. Also relevant is 
the polydispersity p of the simulated sample, defined as 
p = y/ ((di — a) 2 ) /a, where the averages are over all par- 
ticles i. Before any estimation is applied, the best guess 
for each particle size is a with a fractional uncertainty p. 
If the mean estimation error Aa is less than p, the es- 
timation method improves our knowledge of the particle 
sizes; we will show this is true for the simulated data. 

Our first test case is a random close packed sample. 
In such a sample particles do not move, and so time- 
averaging cannot be used. However, particles are packed 
so that they contact each other, that is, 6 = 0. The 
number of contacting neighbors varies from particle to 
particle, so it is not clear how many neighbors should be 
considered. Accordingly, we plot Aa as a function of Z 
in Fig.[IJ4. We find that Aa is a minimum at Z = 5, and 
is indeed much smaller thanp (0.01 vs. 0.07 in this case). 
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It is possible that while Aa is small, that there are 
systematic errors depending on the real particle size a^. 
To test this, in Fig. [2]A we show the ratio between the 
estimated radius and the given radius a\ ■ /o» as a func- 
tion of dj for a polydispersity p = 0.07 RCP sample. 
The symbols and the error bars correspond to the mean 



and standard deviation of the distribution of a^ lu /oj be- 



,(io) 

tween [a.^ai + 0.01], respectively, a- 10 ' /a.; should be 1 if 

our estimation is perfect and indeed we find af 0) / a, = 
1.000 ± 0.013. The quality of the results is nearly uni- 
form as a function of particle size. To check the validity 
of our method for RCP samples with different polydisper- 
sity, we plot the uncertainty Aa as a function of sample 
polydispersity p in Fig. [T]S. We find Aa w p/6 [34| . 

A colloidal gel shares a similarity to a RCP sam- 
ple (touching particles), and has a significant difference 
(much lower volume fraction) . In a colloidal gel particles 
are stuck to their neighbors and form a large network. 
Often the attractive interactions have a finite range, for 
example with depletion gels [37j [see discussion in Meth- 
ods]. Thus we note that the distribution of ay for gels 
is slightly broader than that for RCP, though the mean 
average of Sij is close to 0. Some time averaging is possi- 
ble, although such samples are frequently noncrgodic or 
at best rearrange quite slowly. 

Likewise the contacting particles make gels similar to 
RCP samples locally. However, the contact number fluc- 
tuates greatly in a colloidal gel, and the number of neigh- 
bors averaged over must vary from particle to particle. 
Rather than being a fixed parameter Z, we have a vary- 
ing number of neighbors Zi used in the average (Eq. [5]). 
To determine Zi, we define the coordination number Ci 
as the number of particles within a distance 2.8a, which 
is the first minimum of the pair correlation function. 
We find the average coordination number c « 13.1 for 
a RCP sample, but this will generally be smaller for a 
gel (l6j . Thus for every particle we estimate the num- 
ber of touching neighbors Zi = 5ci/13 where we round 
Zi to the nearest integer. In general, given the tenuous 
nature of a gel, for many particles Zi is fairly small; also, 
Sij has a broader distribution, and so Aa will be worse 
than the RCP case. However, Aa is improved by time- 
averaging, which also minimizes the uncertainty due to 
particle tracking errors. Fig. [2]B shows the ratio between 
the time-averaged estimated radius and the given radius 
°%^/ a i as a function of the true radius a^ for the colloidal 
gel. We find that a^/a, = 1.000 ± 0.018. Aa = 0.018 
is much smaller than the polydispersity p = 0.07. 

Moving from gels, we next consider a dense suspension 
of purely repulsive (hard-sphere) particles. Here no par- 
ticles are in contact, so Sij has a much broader distribu- 
tion; however, time-averaging is even more powerful. We 
show a\ /di as a function of a% at <fi = 0.51 in Fig. [5]C, 
finding a ( i W) /a l = 1.000 ± 0.014. Yet again Aa = 0.014 
is much smaller than the polydispersity p = 0.070. 

For a dense suspension it is not obvious how many 
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FIG. 2: (Color online) The estimated radius o' 10 ^ /di as a 
function of the true radius di, for four simulated systems: 
(A) random close packing, (B) a colloidal gel at cj> = 0.10, (C) 
a nominally monodisperse suspension at 4> = 0.51, and (D) 
a nominally binary suspension at <j> = 0.51. For (A-C), the 
polydispersities are p = 0.07. For the binary sample, the size 
ratio is 1 : 1.3, the number ratio is 1 : 1, and each species 
has an individual polydispersity of p = 0.04. The error bars 
correspond to the standard deviation of of a' 10 ^/ai between 
Oi and a; + 0.01. The insets show sketches of each system. 



nearest neighbors should be used in the average (Eq. [4} , 
so we plot Aa as a function of Z in Fig. [T]A for two differ- 
ent volume fractions. Aa is minimized at Z = 7 for the 
non-RCP samples (circles and triangles in the figure), so 
we fix our choice Z = 7 for all our <p < 0.6 experimen- 
tal data (discussed below). Figure [IJ4 demonstrates that 
Aa does not depend too sensitively on this choice. How- 
ever, it should be expected that for a more dilute system, 
the importance of caging decreases, and the number of 
neighbors a particle has will fluctuate significantly. For 
fixed polydispersity p = 0.070, we find Aa = 0.023 for 
<j) = 0.45 and Aa = 0.060 for = 0.40. This suggests 
that for <f) < 0.40, the estimation method may not be 
useful without further modifications. 

To check the influence of the sample polydispersity at 
fixed = 0.51, we vary p with results shown in Fig. [T]B 
(triangles). We find Aa 0.005 +p/7, suggesting that 
the estimation is useful for samples with p > 0.01, that 
is, any realistic sample. Aa is nonzero when p — 0, in 
contrast to the RCP case. This is due to the distribution 
of Sij in a dense but non-contacting sample. 

The last case we examine with simulation data is a 
nominally binary sample. We simulate a dense suspen- 
sion composed of particles with a size ratio 1:1.3 (mean 
sizes 0.877 and 1.14) and number ratio 1 : 1. For both 
"small" and "large" particles, there is a polydispersity 
p = 0.04. The results are shown in Fig. and we find 
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FIG. 3: (color online) Particle size distributions. (A) The 
distribution of the true radius (solid lines) and the estimated 
radius a' 10 ' (triangles) for a simulated nominally monodis- 
perse sample with (f> = 0.51 and polydispersity p = 0.07. 
(B) The distribution of the true radius (solid lines) and the 
estimated radius a' 10 - 1 (triangles) for a simulated nominally 
binary sample with <j> = 0.51; see text for more details. (C) 
The radii distributions for the nominally monodisperse ex- 
perimental suspension from ref. 0], for five volume fractions 
as indicated. The solid line is a Gaussian fit to the com- 
bined data, giving p — 0.045. (D) The radii distributions for 
the nominally binary experimental suspension from ref. 
for four volume fractions as indicated. Here the solid line 
is a fit to the sum of two Gaussians. With the size ratio 
(1:1.3) and the polydispersity of each species (p = 0.049 for 
the small species, p = 0.050 for the large species), the two 
sub-distributions have substantial overlap. 



of 0) /<Xi = 1.000±0.024 at = 0.51. (Here we have fixed 
Z — 7.) Again, there is no strong dependence on the 
true particle size a%, and in particular the particles in the 
tails of the distributions are estimated with good accu- 
racy. However, the uncertainty Aa for the binary mixture 
is larger than what is found for the nominally monodis- 
perse distribution. This is consistent with the overall 
polydispersity of the sample being larger, p = 0.14. 

An important use of the estimation technique is to 
measure the particle size distribution of a sample in situ; 
we wish to validate this idea with the simulation data. 
To do this, we compare the estimated radius distribu- 
tion P(af °' ) ) with the true radius distribution P(a,) in 
Fig. [SJAj-B. In both the nominally monodisperse sample 
and the nominally binary sample, the estimated distri- 
bution (symbols) is quite close to the true distribution 
(lines). Our results show that the estimated distribu- 
tion is essentially reproduced by convolving the true dis- 
tribution with a Gaussian of width Aa. For a single- 
species sample with a Gaussian distribution of radii with 
polydispersity p, the estimated polydispersity would be 
Aa 2 . Given that for most situations we have 



shown Aa <C p, our technique will only slightly increase 
the apparent polydispersity of a sample. 

One key difference between simulations and experi- 
ments is the boundary condition. Our simulations have 
periodic boundaries. In an experiment, we can not find 
all nearest neighbors of a particle when the particle is lo- 
cated at the edges of an image. This situation is similar 
to colloidal gels, where the number of nearest neighbors 
varies for each particle, and we adopt the same solution 
used there. For each particle, we average over a number 
of nearest neighbors given by Zi = 7ci/13, where is 
the observed coordination number defined before, and we 
round Zi to the nearest integer. The denominator 13 is 
chosen as the number of neighbors in a close-packed sam- 
ple, and the numerator 7 is from the results of Fig. [TJ4. 

Furthermore, we need one more improvement when we 
apply our method to a nominally binary sample. It usu- 
ally happens that we know the mean radii of each of 
the two species, while the number ratio of two species is 
unknown, which means that a is unknown. In this situa- 
tion, we start with a reasonable guess for a! to be used in 
Eq. [3J Then we compute the particle radii and obtain the 
double peak distribution which depends on our guess a'. 
Both peak radii of the trial estimated radius distribution 
should be shifted by (a' — a) from the known mean radii. 
Thus we subtract (a' — a) to adjust the peak positions 
to the known mean radius of each species and we obtain 
the estimated particle size. 

In an experiment we do not have an alternate means 
to determine each particle size and so cannot directly 
verify our results in the way that the simulations allow. 
However, evidence that our method works is shown in 
Fig. \3jP,D. Here, we analyzed previously published ex- 
perimental data from ref. 0] (nominally monodisperse) 
and ref. [l2| (nominally binary). In each case, data from 
several different volume fractions are shown. The size 
distributions agree well for the different volume fractions 
for both the monodisperse and binary cases. Each differ- 
ent volume fraction was a sample taken from the same 
stock jar and therefore should have the same size distri- 
bution, so this is a confirmation that our method works 
well with experimental data. 

We now demonstrate the utility of our algorithm by 
studying colloidal crystal nucleation. The nucleation of 
crystals in a dense particle suspension depends sensitively 
on polydispersity [241 , [2SL l29l] . We examine data of the 
4> = 0.46 sample from ref. [10j, analyzed at longer times 
to examine the crystallization process that was discarded 
from the analysis in ref. [l(| ■ These particles are slightly 
charged, shifting the freezing point to 4>{ rcczc ~ 0.38 and 
the melting point to </> mo it ~ 0.42 [l4| . In this data, we 
confirm that the crystal nucleus appears at the center of 
our microscopic image: this is homogeneous nucleation, 
not heterogeneous nucleation near the wall. 

At each time step, we calculate the number of ordered 
neig hbors N a for each particle using standard techniques 
[Uim [see Methods]. By convention, a crystalline par- 
ticle has N Q > 8 [3, HH . At each time we compute the 
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number fraction of the sample that is crystallized, X(t). 
Figure 0J4 shows X(t) as a function of both individual 
particle size and time, where darker colors correspond to 
larger values oiX(t). Below t = 3000 s, X < 0.2 for all a, 
and essentially all crystal clusters are below the critical 
size (~ 100 particles) [Hj]. At t = 3000 s, a sufficiently 
large crystalline region appears and begins to grow. X 
increases first for particles with a close to the mean ra- 
dius, and these particles continue to be the subpopulation 
that is the most crystallized at any given time. At longer 
times the particles with a farther from a gradually begin 
to crystallize. 

We next consider an alternate way of thinking about 
the same data. Figure|4]S shows the relationship between 
the sample-averaged X(t) (solid black line), the poly- 
dispersity px for all crystalline particles (blue circles), 
and the polydispcrsity p„x for all non-crystalline parti- 
cles (green squares). X starts to increase at t = 3000 s, 
and those particles that are crystalline at that time 
have px ~ 0.03, smaller than the bulk polydispcrsity 
p = 0.045. As the sample crystallizes we observe that 
both px and p n x increase. The growth of px indicates 
that the crystal, while nucleating in a fairly monodis- 
perse region, can grow by incorporating particles that 
are farther from the mean size. In the final state, the 
local polydispersity of the crystalline particles has nearly 
reached the mean polydispersity p. The growth of p n x 
indicates that those particles that are still outside the 
crystal are more likely to be those with unusual sizes. 

The spatial distribution of particles at the end of the 
experiment is shown in Fig. \$!P,D. Figure 2]C shows the 
locations of the crystalline particles, while D shows the 
locations of the non-crystalline particles. Green particles 
have a,i close to a, while the smallest particles are drawn 
blue and the largest drawn red. The cores of the crystal 
regions are composed of the green particles. 

Next, we examine the beginning of the crystal nucle- 
ation process. While many particles are close to the mean 
size, only a few end up being the nuclcation site. To un- 
derstand which ones nucleate, we now focus on the parti- 
cles close to the mean size: radii 1.175 /im < a,: < 1.185 
/im. Among those particles, we define the nucleus parti- 
cles as those that are crystalline particles at t = 5000 s; 
the remainder are non-nucleus particles. We next define 
the local polydispersity Pi{r) of particle i as 

Pi{r) = {{aj - a l ) 2 )/a l (5) 

where the angle brackets () indicate an average over all 
particles j with centers within a distance r from particle 
i. Figures [5)4, B show space-time plots of the mean value 
of Pi(r) for the nucleus particles (A) and the non- nucleus 
particles (B). In all cases, p%(r) is lower close to the par- 
ticles and increases with increasing r. However, notably 
the contours for low pi are at smaller values of r for the 
nucleus particles (A). After t ~ 3000 s, the region of low 
Pi spreads to large values of r for the nucleus particles 
(A), while little change is seen for the non- nucleus parti- 
cles (B). Simultaneously, we show the temporal change of 
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FIG. 4: Relationship between individual particle size a,i and 
crystallization. (A) The contour plot of the fraction of crys- 
talline particles X as a function of a,i and time t. The contour 
lines are numbered by the value of X. The dashed line shows 
the mean radius a. Particles with radii close to a crystallize 
faster. (B) The fraction of the sample X that is crystalline 
as a function of time (solid black line), along with the mean 
polydispersity of all crystalline particles (blue circles) and all 
non-crystalline particles (green squares). The gray dashed 
line corresponds to the bulk polydispersity p — 0.045. (G) 
The crystalline particles are shown at t — 20000 s. (D) The 
non-crystalline particles are shown at t = 20000 s. For G and 
D, the color indicates at for each particle, where blue corre- 
sponds to smaller a, light green corresponds to a* close to a, 
and red corresponds to larger ai. 



N Q in Fig. [S]C for the nucleus particles (solid black line) 
and the non- nucleus particles (dashed gray line). This 
confirms that the onset of crystallization at t sw 3000 s co- 
incides with the expansion of the low local polydispersity 
region seen in Fig. This is all evidence that crystal 
nuclei are formed from regions where the particles are all 
similar sizes. A reasonable conjecture is that nuclcation 
rates are possibly quite sensitive to how well-mixed the 
sample initially is, in this respect of local polydispcrsity. 

We have developed a general method to estimate the 
particle sizes in a dense particulate samples where the 
particle positions are known. Simulations demonstrate 
the validity of our method. This method can be applied 
to any cases where three-dimensional particle positions 
can be found; while we have focused on colloidal sam- 
ples, granular media are quite similar 0, [39} . We have 
demonstrated the utility of our method by examining ho- 
mogeneous colloidal crystal nucleation. While it has been 
known that nucleation is faster for more monodisperse 
samples, we find this is true on a quite local scale. Nucle- 
ation happens in regions that are locally more monodis- 
perse, and crystal growth is proceeds by preferentially 
incorporating particles close to the mean size. 
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FIG. 5: Relationship between nucleation properties and the 
local polydispersity Pi(r). (A) The contour plot of the mean 
Pi(r) as a function of r and t, averaged over all particles which 
are crystalline at t = 5000 s. The numbers on the contour 
lines represent the value of Pi(r). (B) The contour plot of the 
mean Pi(r) as a function of r and t, averaged over all particles 
which are non-crystalline at t = 5000 s. The particles con- 
sidered in (A) and (B) are only those with radii close to the 
mean radii; see text for details. (O) The number of ordered 
neighbors N as a function of time for those particles plot- 
ted in (A) (nucleus particles, solid line) and (B) (non-nucleus 
particles, dashed line). 



Materials and Methods 

We simulate four particle suspension systems, which 
arc random close packing (RCP), colloidal gel, single 
component suspension, and a binary system. The poly- 
disperse RCP sample is generated using the algorithm of 
ref. [ioj . For the three other cases, we perform three- 
dimensional Monte Carlo simulations with hard spheres. 



Additionally for gels, we wish to model colloid-polymer 
mixtures and so we use the Asakura and Oosawa model 
[3^ , This model leads to a pair interaction between two 
hard colloidal spheres in a solution of ideal polymers as 
U(r) = 00 for r < a ih U(r) = -j?k B Tp p [r 3 - 3(<7 i3 - + 
R G fr + 2(tr« + R G ) 3 } for a l3 <r< a tj + 2R G , U(r) = 
for r > Uij + 2Rq. where a\j = (oi + <jj)/2. a is a 
diameter of particle i, k B is Boltzmann constant, T is 
temperature, p p is the number density of polymers, and 
Rg is the polymer radius of gyration. We fix Rq = O.Ict 
and 4> p = Ait /SRqPp = 0.1 where a is the mean diam- 
eter of the hard spheres. For our single-component and 
two-component hard sphere suspensions, particles inter- 
act via U(r) — 00 for r < <Tij, otherwise U(r) = 0. We 
use 1024 particles with the mean radius a = 1 and vari- 
able polydispersity for all simulations. 

The experimental data come from prior experiments 
[l2|. These experiments used sterically stabilized 
poly(methyl mcthacrylatc) (PMMA) particles and im- 
aged them with confocal microscopy. The particle po- 
sitions were located and tracked using standard particle 
tracking techniques [IJ HH ■ Detailed experimental dis- 
cussions are in the prior references. 

We use previously developed order parameters to look 
for crystalline particles and ordered structure [13, [H, HJ. 
For each particle i, we find its nearest neighbors j and 
identify unit vectors fy pointing to the neighbors. We 
then define a complex order parameter qi m using qi m {i) = 
Y]j—i Yim(*ij) where c, is the number of nearest neigh- 
bors of particle i and Y/ m is a spherical harmonic func- 
tion; we normalize this as qi m = qi m /N where TV is a nor- 
malization factor such that J2 m 9lm{i)Q* m (i) = 1 E3- We 
use 1 = 6. For each particle pair, we compute the complex 
inner product d 6 = Y, m 9lm{i)qi m (j)- Two neighboring 
particles are termed "ordered neighbors" if d,Q exceeds a 
threshold value of 0.5. For each particle, we focus on 
N Q , the number of ordered neighbors it has at a given 
time. N % measures the amount of similarity of structure 
around neighboring particles. N l D =0 corresponds to ran- 
dom structure around particle i, while a large value of 
TV* means that particle i and its neighbor particles have 
similar surroundings [38j . 
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